Ab initio strategy for muon site assignment in wide band gap fluorides 
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We report on an ab initio strategy based on Density Functional Theory to identify the muon sites. 
Two issues must be carefully addressed, muon delocalization about candidate interstitial sites and 
local structural relaxation of the atomic positions due to /^"'"-sample interaction. Here, we report 
on the validation of our strategy on two wide band gap materials, LiF and YF3, where localization 
issues are important because of the interplay between muon localization and lattice relaxation. 



I. INTRODUCTION 



When collected into a high-intensity spin polarized 
beam, muons become a powerful probe for many fields 
of physics and other scientific areasiA"— Many appealing 
characteristics have determined the success of muon spin 
rotation/relaxation spectroscopy (/xSR). Firstly, /.tSR sig- 
nificantly widens the classes of materials that can be 
studied if compared to other spectroscopic techniques 
(like NMR and ESR for example) since it can be applied 
to virtually any specimen by simply implanting muons. 
Secondly, during its lifetime the muon spin interacts with 
magnetic orders of either nuclear or electronic origin, and 
provides information on local magnetic fields on a small 
length scale and - when fluctuating regimes are involved - 
on a large frequency window. Thus, muons are mainly 
used as a microscopic magnetometer to probe both static 
and dynamic magnetic orders. Relevant results have also 
been obtained when modeling the effects of hydrogen-like 
impurities in semiconductors, and their reaction kinetics, 
in order to study quantum diffusion.—^— Moreover, in spa- 
tially inhomogeneous systems /iSR gives valuable com- 
plementary results with respects to neutron diffraction. 
This also happens with NMR since implanted muons 
probe the sample from the interstitial region far from nu- 
clei. Nonetheless, the muon localization process and the 
site assignment remain longstanding problems and may 
represent a serious issue in many /uSR experiments since 
the muon position is often needed to extract quantita- 
tive information from /iSR-data. After implantation the 
positive muons (/i'*') usually stop at high symmetry in- 
terstitial sites of the crystal lattice. In metals, a cloud of 
conduction electrons efficiently screens positive charges, 
so that muons leave nearly unperturbed the positions of 
neighboring atoms^i- It is known instead that, due to the 
formation of chemical bonding between the muon and its 
neighboring atoms, the final site is in off-center intersti- 
tial positions in insulators and semiconductors J^^— 

From the very beginning of the rise of /iSR following 
the availability of the first experimental facilities, a lot of 
work has been devoted to the determination of the muon 
sites. A precise characterization of the muon intersti- 
tial sites was indeed possible from accurate experimental 
studies of the Knight shift, of the level crossing resonance 
(LCR) and by inspecting asymmetry relaxation rates as 



a function of applied fields in selected compoundsJ^"— 
Nonetheless, in a large number of cases the muon position 
and its effect on the hosting system cannot be inferred 
solely by experimental knowledge, and a reliable method 
to obtain the muon site in condensed matter would be 
of great value. To this aim, a variety of theoretical and 
computational approaches were usedJ^i^Sr— Successful 
results have been recently obtained for metallic com- 
pounds where an estimation based on the electrostatic 
potential allowed to identify the muon sitesj22, 

In this work we want to show that, among the many 
possible approaches^^i^i^Sii^ based on a first-principles 
method. Density Functional Theory (DFT), already well 
known for its success in studying electronic structure of 
solids, is further a powerful, accurate and effective tool 
to explore the /i"'"-sample interactions on a selected set 
of experimental acquisitions. Our work highlights that 
in insulators the use of a DFT approach is preferable 
since the stronger interaction with the local environment 
makes muon position identification a non-straightforward 
task. 

Here we present our results for two fiuorine com- 
pounds, namely LiF and YF3. They are very useful test 
cases for our computational investigation. Both materi- 
als' ground state electronic structure is well reproduced 
by DFT. Therefore we expect it to provide an accurate 
value for the electrostatic potential and the local atomic 
structure surrounding the /x"*". Among the insulators, 
LiF and YF3 are two well studied cases where a strong 
/i "'"-lattice interaction develops. This leads to the forma- 
tion of a trimer structure with an ionized ^"'" between 
two F nuclei, known as F-/i"'"-F complexfii The dipolar 
interaction between F and /i spins in F-/i"'"-F produces 
a signature in the /iSR signal. This allows an accurate 
determination of the /z"'"-F distance, which we use to val- 
idate our calculations' results 

The paper is organized as follows. After presenting the 
computational strategy for the identification of /x"*" sites 
in Section |TT] we describe the procedure based on the so- 
lution of the Hamiltonian for the /i"*"-F nuclear spin in- 
teractions allowing the determination of the local atomic 
environment from /iSR spectra in Section Hill In Sections 
IIVI and |V] we discuss the outcome of our calculations on 
LiF and YF3 respectively. Finally we draw a summary 
and conclusion of our work in Section IVIl 



FIG. 1. Isosurfaces of the electrostatic potential in LiF for Vfi{r) = 50meV, lOOmeV and 500 meV in a), b) and c) respectively. 
The isosurface in c) represents the localization volume (see text) for the muon in the bulk electrostatic potential. 



II. COMPUTATIONAL DETAILS 

We use the Generalized Gradient Approximation 
(GGA)^ to the DFT and the pseudopotential based 
plane wave method (PPPW) as implemented in the 
QuantumEspresso package. Li, Y and F are somewhat 
difhcult elements for different reasons. F is a first row 
element with a deep potential, Li and Y have shallow s 
semi-core levels. To insure convergence we used Projector 
Augmented Wave approach^ with explicit treatment of 
the s semi-core as valence and a plane wave basis set up 
to 800 eV. The Brillouin zone integration is not a critical 
issue in those wide band gap insulators, a gaussian smear- 
ing with a 8 X 8 X 8 Monkhorst-Pack mesh for the bulk and 
4 X 4 X 4 for the supercell ensures convergence. Since we 
are interested in the electrostatic potential generated by 
the electron and nuclei our calculations may be in princi- 
ples biased by the pseudization of the potential inside the 
atomic core region. For this reason we double checked our 
muon distribution calculation by comparing the results 
we got for the perfect bulk system with full-potential cal- 
culations using the Augmented Plane Wave plus Local 
Orbitals (APW-l-lo)^'^''^'* method as implemented in the 
Wien2K package.— We found that the muon position is 
not affected by the approximation in the description of 
the electrostatic potential in the interstitial region due to 
the pseudization. 

The muon stopping site search problem can be ap- 
proached as the solution for the motion of a parti- 
cle slowing down in an effective potential given by a 
mean field approximation. Here coherently with previ- 
ous approaches^Si^^ we consider that the potential felt 
by the muon V^t is the sum of the Hartree and nuclei 
terms: 



n(r') 



dr' 



^ r R, 



(1) 



We call this the electrostatic approximation and we dis- 
regard possible electron-muon correlation effects. Such 
a potential has minima Vq in the interstitial positions. 
Since muons are light-mass particles it is not trivial the 
identification of the minima with stopping sites. Indeed 
zero point motion (ZPM) effects may play an important 
role in muon localization. LiF and YF3 are good exam- 
ples in this respect. Because of the ZPM not all of the 



minima are stopping sites. If we have more than one 
minimum inside a primitive cell the muon hops between 
neighboring sites if the barrier to be overtaken is lower 
than the ZPM energy. At the end the muon will stop in 
a minimum surrounded by barriers as high as to make 
further hops impossible, or it will share more than one 
neighboring minima positions. In the latter case the cal- 
culation of the mass center for the ground state will be 
necessary to identify the muon position. Inspection of 
in three dimensional systems to understand muon delo- 
calization can be difficult, we need therefore a criterion 
to define the extension of the wavefunction about a min- 
imum. We use the turning point concept saying that a 
wavefunction spreads over the volume of space, we call 
localization volume, satisfying the condition V^(r) < Eq. 
Eq is the zero point motion energy (in short zero point 
energy, ZPE) as the eigenvalue of the ground state solu- 
tion of the Schrodinger equation for the /i^ particle: 



2mn 



(2) 



As a first approximation, as we did in a recent work^l, 
ZPEs can be computed by modeling each minimum as 
an anisotropic harmonic well 

Vir) = im^ {ujIx^ + o^y + ulz^) + (3) 

with eigenvalues given by 

E{nx,ny,nz) = h [ujxinx + 1/2)+ 

+ujy{ny + 1/2) + uj,{n, + 1/2)] + Vo, 



and a ZPE 



Eq = h{uJx +i^y+ (^z)/'2 + Vq. 



(4) 



This method cannot be used whenever the potential well 
surrounding the minimum has an irregular shape. This 
is indeed the case here as sections IIVI and |V] will show. 
Using this procedure we will end up with disconnected 
localization volumes each one corresponding to a given 
stopping site. It may also happens that a localization 
volume goes across the primitive cell boundary connect- 
ing neighboring cells. This would correspond to muon 
diffusion across the crystal, a possibility that, as we will 
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show, is prevented by the formation of the F-/i+-F com- 
plex. 

If the muon did not modify its environment, that pro- 
cedure would be sufficiently accurate. Instead muons in- 
duce a relaxation of the local neighboring lattice struc- 
ture. Computing the local lattice relaxation we neglect 
the spread of the muon wavefunction due to the ZPM, 
then, in the framework of DFT, the effect of muon trap- 
ping on the surrounding structure can be studied as it 
was the trapping process for an interstitial hydrogen. We 
set up a supercell up from our bulk crystal and insert an 
hydrogen interstitial atom. Since we don't want to re- 
produce the actual dynamics of the implantation process 
for the muons we use an heuristic approach to find the 
stopping site. We place the hydrogen impurity in the 
sites identified as minima by the electrostatic potential 
landscape technique or by insights for experiments. Then 
we let the system to evolve to the ground state allowing 
both electron rearrangement and lattice distortion. The 
final optimized position for the impurity represents the 
refined muon position. Since this approach needs an ed- 
ucated guess on the position of the muon it should not be 
considered as totally alternative to the first. Moreover it 
does not include any effect due to ZPM. 

The muon is represented by the hydrogen pseudo po- 
tential in the PAW formalism. We built our supercell 
up from our bulk structure by doubling the bulk primitive 
cell along each crystal axis direction (2x2x2 supercell) . 
A convergence threshold of 5 meV is set for the total en- 
ergy convergence for structural minimization. Here we 
want to study the localization of a /i+ therefore we make 
use of charged supercells. Since charged supercells can- 
not be treated in the periodic boundary conditions we use 
a neutralizing compensating background approach. The 
accuracy of supercell calculations is limited by the size of 
the simulation. For neutral light impurities 2x2x2 su- 
percells may be enough. Since we deal with charged im- 
purities (muon interstitial) we made a convergence test 
using SIESTA codci^ Comparing the structure and to- 
tal energy of 2x2x2 and 3x3x3 supercells we estimate 
the numerical error on the energy and on the optimized 
distances to be ~ 5 meV and ~ 0.02 A respectively. 



III. SOLUTION OF THE SPIN-HAMILTONIAN 

LiF and YF3, are especially useful as test cases for 
our DFT calculations because a precise verification of 
the muon site is obtained by best fitting the asymmetry 
signal produced by the dipolar interaction between the 
muon and neighboring nuclear moments. In fluorides, 
because of the high nuclear moment of F nuclei (^^F has 
spin 1 = ^/2 and ^ 100% natural abundance) and of the 
high electronegativity of this element, the interaction, 
commonly referred as F-/i^-F, is more pronounced. An 
entangled quantum state develops between the muon and 
the surrounding nuclei and the system may be described 
with the following Hamiltonian 
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:S,.S,-3(S,.r)(S,.f)], (5) 



where r is the vector between spins Si and Sj of either 
the fluorine nuclei or the muon, which have gyromagnetic 
ratios 7^ and 7^. The muon depolarization is given by: 



m.n 



(6) 



where TV is the Hilbert space dimension, \m) and \n) are 
eigenstates of % and huJm,n are the corresponding eigen- 
values, fj^ is the Pauli spin matrix corresponding to the 
quantization direction and h is the Planck constant. In 
a powdered sample with cubic symmetry, the observed 
signal is the results of the weighted average over all di- 
rections, i.e. 



|(m|ac|n)P = i(|(m|a,|n)p+ 

+ \{m\<Jy\n)f + \{m\a.,\n)\' 



(7) 



Since the dipolar interaction is inversely proportional 
to the cube of the inter-nuclear distance, one usually con- 
sider only up to next neighboring atoms in order to make 
the calculation of the muon polarization computationally 
inexpensive within a negligible loss of accuracy. More- 
over, the coupling between F nuclear spins may be often 
disregarded with a limited loss of accuracy even if its in- 
clusion does not lead to an increase of the computational 
load. 

For an axially symmetric F-/i+-F complex, as in the 
case of LiF, when considering only two nearest neighbor- 
ing F atoms, the analytic solution of Eq. [5] for a powder 
averaged depolarization is; 



2 / 13 \/3 

Gp{t) = - ( 3 + cos VSujdt + (1 - cos — - — ujdt 

In 3 + V3 \ 

+ (1 + -j=) COS — - — udt I 

(8) 

where oj^ = f^QjFJtih/ {2r^). As will be shown hereafter, 
Eq. [8] fails to capture the data trend in YF3 , and the 
F-//+-F effect alone is not sufficient to determine the /i+ 
site. 



IV. LiF 

LiF has the NaCl crystal structure, with a four for- 
mula unit conventional cubic cell, containing eight cubic 
cages with vertexes at four Li and four F atoms. As 
shown in Fig. [T] the minima of the electrostatic poten- 
tial in LiF are located approximately at the center of 



FIG. 2. Possible muon sites in YF3 (left and center) and LiF (right). The label A identifies the expected site in both compounds. 
Localization volume surfaces are shown in dark yellow for YF3 and in Fig. [He) for LiF. The electrostatic potential section in 
YF3 and LiF allow a direct comparison the the ZPEs in the two compounds. For sake of clarity, unrelaxed lattice structures 
are shown. 
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FIG. 3. Expected /i asymmetry spectra for optimized muon 
sites and atomic coordinates in powdered LiF. Sites are indi- 
cated as shown in Fig.[T] All calculations include all the neigh- 
boring atoms giving rise to couplings higher than one tenth 
of the maximum coupling constant (see Eq. O the number of 
F atoms considered depends on the fi'^ interstitial site). Po- 
sition A gives the best agreement with the measured data.— 



each cage. Minima inside the cage are five. The one we 
label B in Fig. [T{b) is at the very center of the cage sur- 
rounded by four equivalent minima labeled C placed in 
the direction of neighboring F atoms. The minima be- 
come connected for E > 75 meV forming a sort of tetra- 
hedron shaped structure with centroid in site B. All of 
these positions are incompatible with the experimental 
muon site which is known from literature and was ob- 
tained with the strategy explained in Sec IIIII ^^ Fig. [Ijc) 
shows the /i+ localization volume according to the ZPE 
obtained by the solution of the Schrodinger equation for 
the muon in the bulk electrostatic potential for LiF. We 
see that the fi~^ is quite delocalized with the localization 
volume forming a connected network across the crystal. 
The experimental position, labeled A in Fig. [IJb), is at 
the boundary of the localization volume. 

Muon delocalization inside a strongly polar solid is the 
condition under which we expect to have a strong effect 
of /j,^-sample interaction on the outcome of our calcula- 
tions. Indeed allowing atoms relaxation in the minimum 
energy configuration we obtain some large atomic dis- 



placement from periodic bulk for all of the sites consid- 
ered here. A strong modification of the crystal structure 
is found when the muon is added to the interstitial site 
A. While F nuclei are attracted by the charged impu- 
rity, Li atoms are repelled. The distance between the 
muon and its neighboring F nuclei is 1.15 A in excellent 
agreement with the experimental data. Also next neigh- 
boring F atoms are affected by the fj,~^ and are subject to 
a displacement of 0.04 A. The relaxed atomic positions 
correctly describe the formation of a F-fi'^ bonding. The 
relaxed structure for the muon sitting in site B shows 
a similar behavior: the distance between the muon and 
its neighboring F atoms reduces from 1.76 A to 1.57 A. 
We note anyway, that the F-/i+ distance in this case is 
too large to reproduce the experimental fast decay of the 
/iSR signal. Site C constitutes a local minimum for the 
structural relaxation and the hydrogenoid impurity re- 
mains trapped there. Anyway the effect of ZPM here 
is important and because of the large delocalization, the 
muon gets out of the local minimum and reaches site A 
as a consequence of the gradual atomic positions relax- 
ation. This behavior is moreover energetically favored 
if we look at the total energies for the /.t+-sample sys- 
tem given by our DFT simulations. The total energy 
for site A, B and C are reported in Tab. HI We see that 
the inclusion of relaxation effects allows to recover the 
agreement with the experimental findings: site A has a 
total energy which is 0.89 eV lower with respect to site 
B and is thus confirmed to be the muon stopping site in 
LiF. The formation of the F-/z"'"-F complex has impor- 
tant consequences on delocalization in LiF. Indeed 
lattice relaxation breaks the lattice periodicity, while the 
formation of a bond with F enhances the /i+ localization 
hindering its diffusion across the material in agreement 
with the experimental evidence. 

The results of our calculations are confirmed by com- 
parison with experimental data. The expected depolar- 
izations for sites A, B and C are shown in Fig. [3] It is 
clear that the time dependencies of the muon polarization 
for the three inequivalent sites are very different allowing 
us to discard sites B and C. Only site A is compatible 
with the observed asymmetry spectra, while the other 
two locations for the fi'^ give significantly worse fits (site 
C) and non physical values for the local modification of 
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LiF 


YF3 




A 


B 


C 


A 


B 


C 


F-fj,~^ distance [A] 


1.15 


1.56 


1.01 


0.144 


1.134 


1.144 


E, - Ea [eV] 





0.89 


0.54 





-0.64 


0.36 



TABLE I. Results for the structural optimization with /i in 
the interstitial positions A, B and C (see text and Fig. [2]). 
Site A is always the experimental/predicted site. F-fi^ is the 
distance between the fi'^ and its nearest neighbor F atom(s), 
Ei — Ea is the difference between DFT ground state energies 
of the relaxed structures. 



the bonds length and distances between and F nuclei 
(site B). Fitting the experimental results with r^+_F as 
a free parameter in Eq. [SI we find that the distorted crys- 
tal structure obtained from DFT calculations reproduces 
the experimental F-F distanc&ii with ^ 1% precision. 



In order to find the muon sites' positions DFT calcula- 
tions are more necessary in YF3 than in LiF. Firstly be- 
cause too many inequivalent /x"*" interstitial positions are 
available in the primitive cell, and so experimental data 
alone do not allow an unambiguous site identification by 
the F-/i+-F signal. Secondly, the Coulomb potential for 
the unperturbed bulk crystal shows only one minimum in 
(1/2, 1/2, 0) that yields a depolarization which cannot cap- 
ture the experimental asymmetry spectra. Moreover here 
the depolarization signal is only roughly captured by the 
axial F-^~^-F expectations as shown in Fig. There- 
fore in YF3 the uncertainty in the muon site assignment 
can be removed only with the help of DFT calculations. 

Following the same procedure detailed before, we relax 
the structure with the muon in non-symmetric interstial 
positions. Six possible inequivalent interstitial sites were 
found after structural relaxations starting from random 
interstitial positions. The three most energetically favor- 
able in-equivalent sites (shown in Fig. [5]) are all close to 
the localization volume. They all are characterized by 
a slightly distorted F-/i+-F bond with the muon shifted 
perpendicular to the F-F axis forming, for sites A and 
C, an angle of ~ 144° between the two bonds. For site 
B the angle changes to ^ 160° (in Fig. ^ the unrelaxed 
structures are shown for the sake of clarity). The three 
remaining sites, will not be considered in the rest of the 
manuscript, since they are too far from the Coulomb po- 
tential minimum, have higher ground state energies and 
result in depolarization functions which are incompatible 
with the experimental results. 

The depolarizations arising from the relaxed structures 
of sites A, B and C are compared in Fig. [5] The relaxed 
energies with the relevant parameters obtained from the 
DFT structural relaxation, are reported in Tab. U while 
ZPEs and energies for the final muon's positions in the 
bulk Coulomb potential are given in Tab. [Ill 

In order to identify the muon site, all the above results 





Eo 


Ea 


Eb Ec 


LiF 


0.50 


0.50 


0.07 0.00 


YF3 


0.145 


0.10 


1.08 0.22 



TABLE II. Zero point energies (Eq) for the muon in the 
Coulomb potential minimum and energies for the relaxed 
muon sites in the bulk Coulomb potential. All energies are in 
eV and only significant figures are reported. 
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FIG. 4. Fit of YF3 data from Ref. [iO with the conventional 
F-/i^-F model (Eq. [8| and with the depolarization calculated 
for the DFT predicted site in the fully relaxed structure. The 
parameters of the fit are detailed in the text and reported in 
Tab, [ml 



must be considered. Indeed, after the structural relax- 
ation, site B, which has the lowest energy, does not pro- 
vide a correct description of the depolarization function 
(Fig. [6]). Instead a good description of the experimental 
data is obtained when considering the expected depolar- 
ization from site A, as shown in Fig. [31 Moreover, we 
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FIG. 5. Expected asymmetry spectra for optimized muon 
sites and atomic coordinates in powdered YF3. Sites are la- 
beled as in Fig. [21 All calculations include all the neighboring 
atoms giving rise to couplings higher than one tenth of the 
maximum coupling constant (see Eq. [51 . Position A gives the 
best agreement with the measured data. 
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Conv. F-/1+-F 


Site A 


Site B 




0.202(1) 


0.184(1) 


0.188(1) 


Pi 


0.77(1) 


0.76(1) 


0.75(1) 




-0.028(1) 


-0.012(1) 


-0.014(1) 




0.18(1) /is"^ 


0.19(1) 


0.15(1) 


/3 


1.27(6) 


1.45(1) 


1.23(1) 




1.23(1) A 


1.17(1) A 


1.22(1) A 


a 


0.73(1) ^s-^ 


0.97(1) ^ls'^ 


1.00(2) /is" 


xl 


4.7 


2.3 


4.1 



TABLE III. Parameters for Eq. [9] obtained from the best-fit 
to the data of Fig. |4] In the first column the resuhs obtained 
with Gfhf defined in Eq. [8] (already obtained by the authors 
of Ref. 1401) are reported. In the second and third columns 
GfhF is calculated from DFT results (see text), r^-^ is 
the distance between the muon and the first neighboring F 
atom(s) for a given jj,'^ site. Small discrepancies between the 
experimental and calculated tf-^ (that may possibly arise 
from the reduced but non- vanishing ZPM neglected in Eq. [5]) 
are accounted by the SljJ parameter (see text). The final F-fi 
distance is obtained by conveniently scaling all the distances 
between the and the atoms included in the sum of Eq. [5] 
All the scaling factors are smaller than 5%. 
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FIG. 6. Fit of YFs data from Ref. \^ with the conventional 
F-fi^-F model (Eq.[H| and with the depolarization calculated 
from DFT results obtained for site A and B. 



point out that kinetics should favor trapping in site A, 
as this is the site with the lowest electrostatic potential in 
the unperturbed bulk structure and the only one whose 
position is inside the /i^ localization volume of the 
global minimum. Both site B and C have energies higher 
than the muons's ZPE. Thus the probability of finding 
the muon in site B is lower with respect to sites C and A 
and therefore the formation of a F-/i''"-F complex in A is 
more likely. This findings led us to conclude that the F- 
/z"'"-F complex in site A is the maximally populated muon 
site in YF3. As for LiF, the localization region shrinks 
as a consequence of the formation of the bond, allowing 
us to neglect the ZPM when solving Eq. [5] 

The experimental data were fitted according to the 



equation: 

A{t) = Ao[plGF^Fit,SuJ)cxp{-{XF^,Ftf) 

+ P2 exp (-(crt)^) ] + Acalbg 

where Aq is the total asymmetry arising from the sample 
and the sample holder, pi measures the fraction of muons 
reaching the F-/Lt+-F site, p2 = 1 — pi and a account 
for the depolarization in the presence of weak nuclear 
coupling and AcaWg is added in order to compensate for 
the background and for the uncertain calibration of the 
non precessing component. Gf^iF is obtained by solving 
Eq. [S] with the lattice structure obtained from DFT cal- 
culations and SbJ is a parameter that accounts for small 
discrepancies between F-/i"^ calculated and experimental 
distances. 

The parameters obtained from the best fits shown in 
Figs. m and [S] are reported in Tab. lIIII We finally add that 
site C is also compatible with the experimental data and 
therefore we cannot rule out the possibility of a partial 
occupation of this site. 

VI. CONCLUSIONS 

In this work we show that, in materials where a strong 
/i+-system interaction is present, the correct interpreta- 
tion of /iSR experiments requires a combination of exper- 
imental and theoretical investigation. The latter is best 
done by an ah initio approach within DFT. Indeed DFT 
is able to provide on the same footing the electrostatic 
potential we use to find candidate sites for fX^ localiza- 
tion and, by the impurity approach, allows us to get the 
refined atomic structure we use to interpret the exper- 
imental data. We tested our procedure in LiF, where 
we showed that the bulk electrostatic potential fails to 
correctly predict the actual /i'^ site, while upon struc- 
tural refinement we were able to reproduce the formation 
of the F-/z'*"-F complex and its structural details (F-/i+ 
distance). We than extended our investigation on YF3 
where the presence of several candidate interstitial site 
makes impossible the identification of the ji^ position 
from experimental knowledge alone. Comparing the ex- 
perimental data with the refined structure obtained by 
DFT investigation we were able to predict the correct 
location and shape for the F-/i"'"-F complex in YF3. We 
point out that such an approach, we testes on materials 
were /i "'"-system interaction is quite large, is of general 
validity and can be applied on a wide choice of different 
material other than wide band gap insulatorsi^ 
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